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ABSTRACT 



Aims. We report on simultaneous observations and modeling of mid-infrared (MIR), near-infrared (NIR) and submillimeter (sub-mm) 
emission of the source Sgr A* associated with the supermassive black hole at the center of our Galaxy. Our goal was to monitor the 
activity of Sgr A* at different wavelengths in order to constrain the emitting processes and gain insight into the nature of the close 
environment of Sgr A*. 

Methods. We used the MIR instrument VISIR in the BURST imaging mode, the adaptive optics assisted NIR camera NACO, and the 
sub-mm antenna APEX to monitor Sgr A* over several nights in July 2007. 

Results. The observations reveal remarkable variability in the NIR and sub-mm during the five nights of observation. No source 
was detected in the MIR but we derived the lowest upper limit for a flare at 8.59/um (22.4 mjy with A 8 59; „„ = 1.6 ± 0.5). This 
observational constraint makes us discard the observed NIR emission as coming from a thermal component emitting at sub-mm 
frequencies. Moreover, comparison of the sub-mm and NIR variability shows that the highest NIR fluxes (flares) are coincident with 
the lowest sub-mm levels of our five-night campaign involving three flares. We explain this behavior by a loss of electrons to the 
system and/or by a decrease in the magnetic field, as might conceivably occur in scenarios involving fast outflows and/or magnetic 
reconnection. 

Key words, black hole physics - Galaxy : center- Galaxy : nucleus- accretion, accretion disks 



1. Introduction 

Due to its proximity compared to other galaxies, the center of 
our Galaxy is the best place to study the nucleus of a spiral 
galaxy and its strong influence on its neighborhood. In the past 
two decades, the advent of NIR adaptive optics at 10m class tele- 
scopes has made it possible to track the movements of stars in 
the innermost regions of the Galaxy. The trajectories of these 
stars, which are perfect Keplerian ellipses with a common focus 
to measurement precision, have demonstrate d there is a super- 
massive black hole at the center of the Galaxy (lEckart & Genzel , 



massive black holes in galactic nuclei, making it the best labo- 
ratory for testing the supermassive black-hole paradigm and for 
studying its unique astrophysical environment. 

One of the peculiarities of S gr A* is its variabi lity at all 
wavelengths (see a review in e.g. iGenzel et al 1 120101) . Distinct 
outbursting events in the lightcurves, called flares, have been de- 
tected on different timesca les from a few minutes to a fe w hours, 
observed in X- rays (e.g iBaganoff et all 1200 ll l2003h and the 
NIR (e.g., IGenzel et all 120031: IGhez et all 12004 IClenet et all 



2003) 



1996L 1 1 997t lOhez et all 1 19981 12003L I2005L ISchodel et alTT2002 



At a distance of ab out 8.33 + . 35 kpc, this object has a mass 
of about 4 x 10 6 M© (IGhez et all [20081: iGillessen et all 120091) 
located at the positio n of the compact radio source named Sgr A* 
dBower et al.U2004t) . The Schwarzschild radius of the black hole 
has an angular size of about 10 //as as viewed from Earth. This 
is three times large r than the black hole s in M87 (3.66 Atas , 
iBettoni et al.l l2003h and in M31 (2.28 ^as. iBender elatl 120051) . 
The GC represents the best evidence so far that there are super- 



Send offprint requests to: X. Haubois 

* Based on multi-wavelength observations obtained at the ESO VLT 
Melipal-Yepun telescopes and at the APEX facility (run ID: 179.B- 
0261). 



120051) . These X-ray/NIR flares are also sometimes followed by 
a sub-mm flare with a time delay of 100 min to a few hour s 
(lEckart et all 12001 l^isef-Zadeh et all 120091: iTrap et alll20Tlh . 
Flares ar e strongly polarized at the highest flux levels (for a re- 
view, see lGenzel et all 1201 Of) . Thank s to their detailed polariza- 
tion analysis of some Sgr A* flares, IZamaninasab et all (1201 Ol) 
support the idea of a compact source for the emission, instead of 
radially extended shapes. 

Moreover, these flares display a quasi-periodicity of about 
20 minutes in the NIR which has been interpreted in terms of 
the orbiting motion of a hot spot around a rotating (Kerr) su- 
permassive black hole (IGenzel et all 20031: iTrippe et ail, 120071 : 

l2006ah r 



iMever et all 120061: lEckart et al 



of Lorentz factor y ~ 10 j - 10 6 are required to produce 
the NIR and X-ray emission via synchrotr o n emi s sion and/or 
inverse Compton scattering dMarkoff et all 120011 : lYuan et ail 
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Fig. 1. Chronogram of the observations 



The lightcu rves were created with ape rture photometry in the 
same way as in lDodds-Eden et all (1201 lh . using two small aper- 
tures, one centered on Sgr A*, the other centered on the close 
contaminating star S17. At the time of these observations, Sgr 
A* was conf used with the star S17, a dding a 3.8 mJy to the 
Ks-band flux dDodds-Eden et all 1201 lb . but only a contribution 
of ss 1 mJy to the L'-band flux. We subtracted the contribution 
due to S17 from the lightcurves shown in Fig. [3] However, we 
note that there may be an additional contribution from faint stars 
dDodds-Eden et all 1201 lh . and in L'-band, perhaps some small 
contrib ution from the dust cloud Sgr A*-f that is very close to 
Sgr A* (IClenet et alll2005b . We concentrate only on lightcurves 
in this paper and leave the polarimetric properties for a future 
paper. 

The infrared extinct ion towards the G C has recently been 
investigated in detail by iFritz et all (1201 lh . who derived an ac- 
curate extinction law from hydrogen lines, calibrated with a 
(extinction-free) radio map. We thus used the IFritz et all (1201 lh 
updated extinctions of Au = 1.09 and Ag s = 2.42 for a red 
(compared to a stellar source) Sgr A*-like source. 

2.2. MIR 

MIR observations of Sgr A* at a wavelength of A = 8.59 /mi 
were carried out on 20, 22 and 23 July 2007, using the instru- 
ment VISIR dLagage et all |2004|) . on UT3 at the VLT. These 
observations were taken in parallel with the NIR observations 
in a "trigger" mode; i.e., whenever increased NIR emission was 
seen in real time with NACO, VISIR began immediately to ob- 
serve Sgr A* to increase the likelihood that a flare would be 
detected in both spectral bands. On July 22 for example, our first 
VISIR data point on Sgr A* was obtained around 0:23 UT, about 
30 min after the first (L'-band) peak, and at the beginning of the 
brightest, second peak. 

We observed in the PAH1 filter centered on 8.59 /mi, the fil- 
ter with the highest transmission available on VISIR according 
to the user manual Q. At these wavelengths, the thermal back- 
ground emission dominates the images. We used the classical 
chopping/nodding technique with chopping and nodding throw 
angles of 10". We chose this parameter because we wanted to 
continuously monitor the GC on the detector in order to increase 
the signal-to-noise ratio and to avoid gaps in the temporal cov- 
erage. The principle of this chopping/nodding technique is to 
subtract the high-level sky background (by chopping the sec- 
ondary mirror) and correct the differential instrumental back- 
ground caused by the chopping (by nodding the telescope, de- 
tails can be found in the VISIR user m anual). 

VISIR was used in BURST mode dDoucet et al] . l2006h . The 
principle of the BURST mode is to take shorter exposures than 
the atmospheric coherence time. The main asset of this technique 
is its spatial resolution. In principle, it allows one to reach the 



I2003L 120041: Ibodds-Eden etafl l2009t lEckart et all l2006bh . The 
possible origins of the hot electrons include magnetic reconnec- 
tionin the accretion flow dBaganoff et all 120011 : iMarkoff et all 
1200 lh . stochastic acceleration dLiu et all 120061) . fluctuations in 
the accretion rate, or even the tidal disruption of a small infalling 
body such as an asteroid dCadez et all 120081) . 

While many papers conceive of flares as independent events 
that are limited in time, other teams claim that the IR emission 
from Sgr A* is made up of permanently fluctuating emission that 
could be described by red noise statistics. The p roponent s of thi s 
picture conside r the term "flare" as impr oper (iDo et all l2009h . 
More recently, iDodds-Eden et al ] (1201 lh has resolved the dis- 
crepancy in these two viewpoints, finding that the NIR emission 
has a dual nature, consisting of two components: the first a low- 
level, permanently fluctuating component, identified as the NIR 
"quiescent" state, the second consisting of additional, sporadic, 
high flux "flare" events. 

Revealing the enigmatic origin of the flares of Sgr A* is 
one of the aims of an ESO large program (LP) of observa- 
tions (Proposal ID: 179.B-0261, PI : Pr. R. Genzel). This project 
combined all VLT instruments available for high-resolution in- 
frared studies, as well as sub-mm observations with the Atacama 
Pathfinder Experiment (APEX), and was part of a large ground- 
based and space-based, international multi-wavelength cam- 
paign. Constraining the SED of flares by photometry at a high 
sampling rate in sub-mm, MIR, and NIR can constrain the emis- 
sion process responsible for the flares. A part of this LP was 
therefore devoted to observing Sgr A* simultaneously in these 
three wavelength ranges. Throughout this paper we use /? for the 
spectral index that is defined through vL v oc \fi . For synchrotron 
emission this can be related to the particle index, p, of the 
synchrotron-emitting electron distribution defined as n{y) oc y~P 
through/? = (3 -p)/2. 

In this paper, we analyze simultaneous observations of Sgr 
A* at NIR, MIR, and sub-mm wavelengths, carried out over five 
consecutive nights in July 2007. We begin with the observations 
and data reduction in Sect. [2] In Sect. [3] we present our results. 
After discussing the implications on the spectral energy distribu- 
tion (SED) of our observations in Sect.|4] we conclude in Sect. [5] 



2. Observations and data reduction 

In this section we present simultaneous multi-wavelength ob- 
servations carried out during July 2007 at NIR (2.18 //m and 
3.80 /mi), MIR (8.59 /mi) and sub-mm (850 /mi). A schematic 
log of the observations is presented Fig. Q] 

The July 2007 campaign ran for five nights from 18 July 
2007 to 23 July 2007. On most nights NACO operated in polari- 
metric imaging mode in /iT-band, with some short observations at 
other wavelengths (H, L'-band), except for July 22 when most of 
the night was observed in L'-band. A computer crash at APEX 
on July 18 stopped the sub-mm observations short at ~2:00 UT. 
On the night of July 21, technical problems with NACO pre- 
vented us from obtaining any data, so this night is not included. 

2.1. NIR 

Here we describe the Ks and L'-band imaging ob servations ob- 
tained with NAOS-CONIC A on UT4 at the VLT dLenzen et all 
I2003L Rou sset et all I2003h . The data reduction was the same as 
in lDodds-Eden et all d2009h when we subjected the raw data to a 
sky subtraction computed from jittered object images, followed 
by flat-fielding and a correction for dead/hot pixels. 



1 http://www.eso.org/sci/facilities/paranal/instruments/visir/doc/index.html 
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diffraction limit in good observing conditions (0.22" at 8.59 mi- 
crons, about three times the pixel scale, with an optical seeing 
<0.5"), since the atmospheric turbulence is frozen at the rate at 
which the images are recorded (an integration time of 20 ms in 
our case). This gain in angular resolution increases the sensitiv- 
ity since we have a sharper point spread function (PSF) than with 
the usual mode. 

In BURST mode, a typical data set is composed of two 
nodding cubes (nodding position A and nodding position B). 
One cube contains 3000 images of the two chopping positions 
recorded by sequences of 100 frames. The exposure time of 
one single frame is 20 milliseconds. Thus, a whole cube spans 
one minute in total. For reduction, the raw frames were first di- 
vided by a flat-field calibration image. Using the classical chop- 
ping/nodding correction, one obtains four quadrant images with 
the target appearing as a positive source in two quadrants and a 
negative source in the other two. However, due to variations in 
the atmospheric conditions, the image quality varies from frame 
to frame and a fraction of them are unusable. We thus selected 
better frames, based on the sharpness of the PSF as estimated 
on IRS21 (which was the only available star to do so in the field 
of view we observed). Altogether, 85% of the data are combined 
in a single cube to increase the signal-to-noise ratio. Photometric 
calibration was performed using the observations of several stan- 
dard stars obtained with our scientific data and the stellar fluxes 
from lCohen et al. ( 1999). The final images correspond to about 
three minutes of accumulated integration time, and they have 
been co-added in Fig.f2]to show the detected structures better. 

We performed aperture photometry at the position of Sgr A* 
on all the images taken during the three VISIR half nights to 
search for variability. The position of Sgr A* was determined 
using the pos itions of the GC IRS stars listed in iPaumard et alj 
(120061) and in lOttet al.l (Il999l) . We applied a low-frequency fil- 
ter to these lightcurves to detect variability more easily with time 
scales of about five minutes. However, the use of the classical 
chopping/nodding data reduction was unsuitable for constrain- 
ing the lightcurve of Sgr A*. Indeed, the lightcurves obtained 
are different for each quadrant and display a highly noisy pat- 
tern because of contamination by field sources (faint stellar or 
extended diffuse objects) and atmospheric turbulence. As a re- 
sult we did not perform the classical nodding correction, esti- 
mating instead the differential instrument emission (introduced 
during the chopping phase) from background zones and subtract- 
ing it from the Sgr A* signal. In our case, this technique gives a 
(sky+instrument) background-corrected flux estimation, free of 
bias from contaminating sources. 

For the 8.59^ m extinction we used Ag, ^ Mm = 1.6 + 0.5, from 
iFritz et all (1201 ll) computed for this specific VISIR filter and a 
red (compared to stars) source spectrum. 

2.3. Sub-mm 

The 870/itn data were t aken with th e Large APEX bolome- 
ter camera, LABOCA (ISiringo et al.l 120091) . located on the 
Atacama Pathfinder Experiment (APEX) telescope at Llano 
Chajnantor, Chile. LABOCA is an array of 295 composite 
bolometers operating in the 345 GHz transmission window. Its 
bandwidth is ~ 60 GHz and the FWHM of the PSF at 870/rni is 
19". The atmospheric precipitable water vapor content (PWV) 
during the observations was between 0.6- 1.2 mm corresponding 
to a zenith opacity of 0.3 - 0.45 for the LABOCA passband. 

The Galacti c Center was observe d using the raster-spiral ob- 
serving pattern (ISiringo et al.l 120091) . In this mode the telescope 



performs a spiral at constant angular speed at four positions 
(the raster positions), leading to a fully sampled map over the 
full 11 ' field of view of LABOCA. The telescope scanning 
speed in this mode is between 1'/* and 2.5' /s. Each spiral takes 
about 35 s yielding a total integration time of 140 s and a typical 
rms noise level of ~ 60 mJy/beam for each map. The Galactic 
center observations were followed by observations of G10.62, 
a standard LABOCA secondary calibrator, and G5.89 0. The 
atmospheric zenith opacity was determined via skydips every 
one to two hours. 

The data were reduced using the BoA software package 0. 
Reduction steps on the bolometer time series include temper- 
ature drift correction based on two "blind" bolometers (whose 
horns were sealed to block the sky signal), flat fielding, calibra- 
tion, opacity correction and correlated noise removal on the full 
array as well as on groups of bolometers related by the wiring 
and in the electronics, flagging of bad bolometers and despiking. 
Each reduced scan was then gridded into a spatial intensity and 
weighting map. These reductions step were applied to the Sgr A* 
and the G5.89 and G10.62 scans. Calibration was achieved using 
G10.62 which has a flux density of 33.4 + 2.0 Jy for LABOCA 
dSiringo et aUl2009l) . The absolute calibration accuracy is about 
10%. To estimate the relative calibration error we applied the 
calibration curve of G10.62 to G5.89 and determined the disper- 
sion of the measured G5.89 fluxes over the observing period of 
7.5 hours. This yields a relative calibration accuracy of ~ 2%. 
For the determination of the Sgr A* light curve we first gener- 
ated a model of the sub-mm emission in the GC region by co- 
adding all calibrated, pointing-drift-corrected GC maps. From 
this high signal-to-noise map the point source Sgr A* was fitted 
by a Gaussian and subtracted from the co-added map. We then 
reduced all GC scans again and subtracted the model signal from 
the time series of each bolometer. The resulting maps only con- 
tain the point source Sgr A*. The lightcurve was constructed by 
fitting a Gaussian to each map. To test the relative calibration 
accuracy we also measured the residual peak flux at a position 
near Sgr A* located on the GC mini-spiral. The residual varia- 
tions were found to be 1.5% compared to the model flux at this 
position. The LABOCA light curve is shown in Fig. [3] 

3. Results 

3.1. High-quality map of GC region at 8.59fim 

The first interesting outcome of these data are the very high- 
quality images of the GC area at 8.59 microns. In Fig. [2] we 
see the MIR emission from both the photosphere and the sur- 
rounding dust shell of the bright stars GCIRS21 and GCIRS 1W. 
Also prominent at MIR wavelengths are the gaseous and dusty 
structures such as the Northern and Eastern Arms of the mini- 
spiral. To see both Sgr A* neighborhood and bright stellar and 
gas components, we present the field of view in two different 
images. One of the most striking features in these images is 
the presence of elongated filaments in the Northern Arm. A 
proper moti on study of t hese s tructures in L'-band (3.8 fim) is 
reported in iMuzic et aH (I2007I) . But even more striking is the 
presence of a stretched emission zone very near Sgr A* in pro- 
jection. This structure seems to begin from the Northern Arm 
(2" southeast of Sgr A*) and follows an asterism consisting of 
a few aligned stars leading to the S-star region (those stars are 

2 See http://www.apex-telescope.org/bolometer/laboca/calibration/ 

3 BoA: http://www.astro.uni-bonn.de/boawiki/Boa 
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Fig. 2. Image of the GC zone taken in the PAH1 filter resulting from the co-added images of 3 hours of observation on 22 July 2007. 
Upper panel: Sgr A* surroundings compound of stellar sources and gas filaments. Lower panel: close-up on the Sgr A* zone. No 
extinction correction was performed at this point. 



detected in NIR bands). Then, in the most central part (at around 
0.5" from Sgr A*, see lower panel of Figf2), it deviates a bit 
to the west, then turns back northwest at the Sgr A* position 
(Fig.0. This last segment of the structure is also reported in pre- 
vious papers as a dust ridge whose emission could contaminate a 
potential poin t -like e mission from Sgr A* (IStolovv et al.L Il996t 
ISchodel et all I2007I) . This elongated structure might be due to 
gaseous material falling towards Sgr A* from both sides (south- 
east and northwest), a bridge of material in projection along the 
line of sight or a combination of foreground/background compo- 
nents that appear as a unique element in projecti on. We note that 
the L'-band source Sgr A*-f (IClenet et alll2005l) could belong to 
this filament, but we cannot conclude anything about the nature 
of the structure here. 



3.2. Variability of Sgr A* 

The NIR and sub-mm lightcurves are shown in Fig. [3] In this 
section we discuss the variability in the NIR (Section [3.2.U , an 
upper limit for detection in the MIR (Section r3.2.21 >, and then the 
sub-mm variability and its relationship with the NIR variability 
for this campaign (Section[3723]l. 



3.2.1. NIR variability 



The NIR lightcurves (Fig. [3]l display significant variability on 
all five observation n ights. If transla t ed to Ks-band using a 
color of F v oc y ° 6 dHornstein et al.L l2007l) . the L'-band ob- 
servations would be fainter by a factor «0.7; the //-band ob- 
servations would be a factor =s 1.2 brighter. For compari- 
son, we computed some statistical quantities to describe the 
variability of the lightcurves in Table Q] The strongest flare 
is seen on July 20 (ZTs-band; ^7 mJy). On July 22 an in- 
teresting variability pattern is seen in the L'-band, with four 
clear peaks separated by ~40 min. This fla re was also inves- 
tigated previously by lHamaus et al.l (120091) . who derived an 
average period of 45 minutes. This is more than twice as 
long as th e subpeak-to-subpeak timescales observed in previ- 
ous flares (iGenzel et all 120031: iMever et all 120061: lEckart et all 



l2006bl:lTrippe et all 120071 Dodds-Eden et al 1 120091)" Moreover, 
we point out that the flux drops between subpeaks are also very 
deep. Another flare is seen in the /fs-band just after 5:00 UT on 
July 23. 
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Fig. 3. Lightcurves of Sgr A* at NIR and sub-mm wavelengths for five days of the July 2007 campaign. The green, red, and 
blue datapoints of the upper panels represent the H, Ks, and L'-band NIR fluxes, respectively. Light red squares show Ks-band 
observations for which max(FWHM v , FWHM V ) > 15 pixels (^200 mas). The sub-mm lightcurve is shown below in black. One 
should bear in mind that the NIR observations were not all made in the same filter. 



Table 1. Statistics of the NIR lightcurves. 



NEAR-INFRARED SUB MM 



Night (July 2007) 


18 


19 


20 


22 


23 


18 


19 


20 


22 


23 


Median + Median deviation (mJy, Jy) 




1.41+0.50 






2.96 + 0.14 




Median (mJy, Jy) 


0.99 


1.15 


1.38 


1.52 


1.35 


3.09 


2.96 


2.75 


3.10 


2.91 


Median deviation (mJy, Jy) 


0.27 


0.61 


0.46 


0.56 


0.67 


0.09 


0.11 


0.15 


0.11 


0.18 


Duration of the sample (min) 


338 


359 


351 


344 


500 


174 


422 


388 


444 


506 


# peaks above 


1 


1 


2 


2 


1 

















3-cr threshold of entire obs* 






















# dips below 























2 








3-cr threshold of entire obs* 






















Maximum (mjyjy)* avg of the 3 highest points 


3.21 


3.85 


6.72 


5.79 


5.21 


3.43 


3.25 


2.98 


3.40 


3.32 


Minimum (mjyjy)* avg of the 3 lowest points 












2.99 


2.71 


2.27 


2.79 


2.40 


Average timespan 






120 


45 


40 






120 






Between peaks/dips (min) 























Notes: For all nights, the NIR observation were done in the K-bmd except for July 22 where the L' filter was used. We do not consider the 
/f-band data here. * A peak/dip is defined as a set of more than 5 points in a row (3 points for sub-mm data) above/below the 3-cr level. Flux 
units are in mJy for NIR, Jy for sub-mm. 



3.2.2. MIR upper limit 



The difficulty of detecting Sgr A* in the MIR (compared to the 
NIR) comes from to the lower spatial resolution and the dom- 
inating background. Despite the trigger mode used to observe 
with VISIR only when NIR flares were detected, we did not de- 
tect any point-like source at the position of Sgr A* at 8.59 fim 
in our observations, nor was there any significant correlation be- 
tween the NIR and MIR lightcurves obtained via aperture pho- 
tometry. Up to now, Sg r A* has never been detected in MIR and 
only upper limits exi st ( Telesco et allll995l : IStolovv et allll996t 



ISchodel et all 120071 1201 ]]). The lowest upper limit at 8.59//m to 
date in the literature gives a mean flux estimation of the Sgr A* 
emission of 45+/-13 mJy (ISchodel et all 1201 ll) . Once dered- 
dened with t he upd ated 8.59yum extinction of Ag 59^,,, = 1 .6 ± 0.5 
dFritz et all 1201 ll) computed for the VISIR PAH1 filter, this 
value becomes 31+/-9 mJy and thus 58 mJy at the 3-sigma level. 



We determined an upper limit to the flux by injecting a 3x3 
pixels Gaussian-shaped source of known intensity into our MIR 
images. For each VISIR observation night we subtracted col- 
lapsed images taken during intervals of low NIR emission from 
collapsed image from the flaring state intervals (with the artifi- 
cially source injected). The detection limit is then derived when 
this artificial point source is detected at a 3-cr significance. 



Because of the high quality of our observations, we are able 
to derive an especially tight upper limit on a transient source at 
the position of Sgr A* of 22.4 mJy (dereddened, for July 22), 
simultaneous to the L'-band flare. This is the lowest upper limit 
ever deriv ed for a flare a t 8.59// m. We emphasize here that, com- 
pared to ( Sch odel et all 1201 ll) who also used this dataset, the 
upper limit reported here is due to the specificity of our method, 
which consists in estimating the lowest detectable MIR variabil- 
ity during an NIR flare. This result thus represents an upper limit 
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on the MIR flux relatively to the NIR flux during a flare whereas 
the Schodel et al. one is an absolute MIR upper limit. 

Combined with the NIR mean flux in the same time interval, 
we could also put a lower limit on the spectral index between 
MIR and the NIR bands. If f3 is defined through vL v ~ v 8 , the 
best constraint is obtained for the night of July 23 with a lower 
limit of/3 > -0.5 (Table©. 

This constrain t is no t as strong as (for example) derived by 
iDodds-Eden et al.l ( 120091) . who derived a lower limit on the spec- 
tral index of fS > 0.0 from simultaneous NACO (3.8/^m) and 
VISIR (11.88yum) observations. The weaker constraint on the 
spectral index found here, despite the high sensitivity of the MIR 
upper limit in this paper, is partly due to the relative faintness 
of the flares ob served during this observ ation campaign com- 
pared to that of IDodds-Eden et ail (12009b and partly due to the 
updated extinction law. In fact, the updat e d exti nction also af- 
fects the lower limit of IDodds-Eden et al 1 (120091 such that the 
constraint should be updated to (3 > -0.3. Still, both lower limit 
measurements leave very little room for a negative vL v MIR-NIR 
spectral index. 

A cutoff at high energy in the flare spectral energy distribu- 
tion could be produced if the acceleration process is inefficient at 
accelerating electrons to these energies. If the acceleration pro- 
cess shuts off early, then it can also be due to electron losses from 
synchrotron cooling. Our measurements clearly rule out that the 
three flares arise from the high-frequency emission of a spec- 
trum from an electron distribution with an exponential cutoff at 
sub-mm-emitting energies (i.e. the thermal electrons known to 
produce the sub-mm emission). While it is hardly disputed that 
the NIR emission {A < 4p.m) from Sgr A* requires nonthermal 
electrons, it is worth illustrating the precise restrictions our data 
place on the presence of an exponential cutoff in the electron 
distribution. Assuming such an exponential cutoff in the elec- 
tron spectrum at an energy of y 'cutoff, with particle index below 
the cutoff of po, 



n{y) ~ y 10 X exp(- 



y 



-) 



Jcutoff 



(1) 



the particle index, p = - dl "j" <7 ^ w ju then be p — p$ H — 



leading to a spectral index/? = (3 - p)/2 = 1/2 x (3 -po - - — ). 



The energy of NIR-emitting 
dRvbicki & Lightmanl[T979l) : 



electrons 



7 NIR - \ h'NIR X 



16wc 
3qB ' 



is 



(2) 



Thus our measurement of p^m > -0.5 for the night of the 23rd 
implies, 



1 7nir . „ , 
-(3 - po ) > -0.5, 

i y cutoff 

such that, 

ycutoff = 8 e > 200(^)- 1/2 



(3) 



(4) 



with 6e = kT/mc 2 , for a thermal distribution (p=2 and y cu toff = 
or 



y CMtof f > 1200(^r 1/2 



(5) 



for a more "bottom-heavy" energy distribution below y cutoff 
(p=3). Thus it is clear that our data require electron energies 



up to maximum y factors of at least y cutoff ^ 200, and do not 
allow the NIR flares to arise from a purely thermal distribution 
of electrons peaking in the sub-mm (9 e ~ 10). 

More important, most measurements of the NIR stellar slope 
between 3.8/mi and 1.65/im have been independent of photo- 
metric calibration and extinction, since the spectral index was 
meas ured only relative to stella r colors of the close star S2. 
Both lHornstein etalJ (120071) and lGillessen et al.l d2006l) . for ex- 
ample, found positive slopes of « 0.4 + 0.2 for high fluxes us- 
ing an 'off-state' background subtraction method, which is con- 
sistent wi th the lower limit on the slope reported here, as well 
as that of Dodds-Ed en et al. (2009). At fainter flux es, however, 
lEisenhauer et al.l (120051) and iGillessen et al.l (120061) find a red- 
dening of the spectral index, up to /? « -3 + 1 (depending on 
the background subtraction method) for dim emission of about 
2-3 mJy. Our results are, however, inconsistent with results with 
such steep spectral indices at low fluxes n. For the three methods 
of background subtraction presented in IGillessen et al.l (120061) 
our results only agree with the "off-state" background subtrac- 
tion method for the dim state, perhaps implying that the other 
methods underestimate the spectral index. However, we cannot 
compare the results at low fluxes so easily because our flaring 
intervals include a mix of both bright and faint NIR emission. 

3.2.3. Submillimeter & correlation with NIR 

A striking aspect of the lightcurves presented in Figj3] is the 
marked anticorrelation seen on the night of July 20, where two 
peaks in the NIR emission are coincident with the two deepest 
dips in the sub-mm emission. In particular it is noteworthy that 
the brightest NIR peak of the five nights of observation is coin- 
cident with the lowest sub-mm flux. Similar behavior (an NIR 
peak accompanying a sub-mm dip) is seen on July 19, as well as 
on the night of July 23. 

There are hints that, a peak in the NIR emission is simulta- 
neous with a peak in the sub-mm emission particularly on July 
18. To test whether the apparent anticorrelation of sub-mm and 
NIR emission in this dataset are significant, we binned the NIR 
flux to the sampling of the sub-mm lightcurve and compared the 
bins for non-zero sub-mm and NIR fluxes (Fig. |4). Computing 
the usual (Bravais-Pearson) correlation coefficient of these two 
variables, an anticorrelation with a coefficient of r « -0.5 is 
significant for the A's-band dataset, with a rejection probability 
of 6.0 x 10~ 6 . The anticorrelation is not significant for the L- 
band dataset. The shape of the correlation plot shown in Fig. [4] 
is mostly caused by the anticorrelation of the July 20 flare. If 
we do not include these data, the correlation is no longer sig- 
nificant and raises the possibility that July 20 lightcurves show 
an extraordinary behavior, associated with a flare event and not 
characteristic of the quiescent emission. 

A similar relationship betw een NIR/X-ray and longe r wave- 
length emission was noted by IDodds-Eden et al.l (1201 Oh in ob- 
servations from April 2007, where a dip in millimeter emission 
was observed to accompany an NIR/X-ray flare. In that case the 
brightest NIR/X-ray flux of the campaign (which was indeed a 
very bright flare, brighter than any in our July campaign) was 
followed by a ~200 min dip in the millimeter emission. The 
millimeter flux during the dip was the lowest observed during 



4 Although our observed fl ares of July 20 and July 22 had peak fluxes 
comparable with the flare of IGillessen et aTJ (120061) . the interval over 
which the spectral index was measured had an average NIR flux of only 
0. 1-0.2 the flux of S2 (and would thus correspond to the "dim state" of 
IGillessen et~afll2006l) . 
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Table 2. Upper limits on the detection of a synthetic MIR flare, as well as the corresponding limits on the MIR-NIR spectral slope 
and particle index. 



Night (July 2007) 


20 


22 


23 


On state time interval/s [UT] 


23:07.18 to 00:13:30, 


0:23:12 to 0:39:48, 


5:34:06 to 5:49:0, 




& 00:50:54 to 2:44:54 


& 1:05:06 to 1:24:48, 


& 6:04:30 to 6:19:30 






& 1:57:00 to 2:00:00 




Off state time interval/s [UT] 


00: 17:30 to 00:35.06 


0:41:06 to 1:03:48, 


5:34:06 to 5:49:06, 






& 1:26:12 to 1:55:36, 


& 6:04:30 to 6:19:30 






& 2:01:12 to 2:39:12 




Upper limit for transient emission 


24.2 


22.4 


31.2 


(dereddened with Ag.s^ = 1.6 ± 0.5) 








Average NIR flux (wavelength) 


2.0 (Ks) 


3.4 (L') 


4.1 (Ks) 


on state time intervals 








Lower limit on the vL v spectral index 


P >-0.80 


P>-\3 


P >-0.46 


(vL v oc v 8 ) 








Upper limit on the particle index, p 


p < 4.6 


p < 5.6 


p < 4.0 


(N(y) oc y-P) 









3.4 
3.2 

I 3 

a 2.8 



M 2.6 
3 



2.4 



2.2 




01234567 
NIR flux [mjy] 

Fig. 4. Comparison of simultaneous NIR and sub-mm fluxes for 
the July 2007 campaign. The graph shows the sub-mm vs. Ks- 
band (red squares) or L-band (blue crosses) NIR flux. The July 
20 night data are marked as filled-in red squares 

. The NIR fluxes were binned to the sampling of the sub-mm 
lightcurve (5.6 min). The Ks-band scatterplots shows a significant 
anticorrelation of r x -0.5, with p-value p = 6.0 x 10~ 6 . 



that entire ~12 day cam paign. The explanation put forward in 
iDodds-Eden et al.1 d2010h was that a dip in millimeter emission 
may occur due to the loss in emissivity of mm-emitting electrons 
that would occur in a magnetic reconnection event, in which sig- 
nificant mag netic energy i s lost in the inner regions of the ac- 
cretion flow. lEckart et al.1 (120081) report a bright NIR flare and 
a sub-mm signal in which the lowest sub-mm flux is also co- 
incident with the highes t NIR flux (between 5h and 6h UT). 
lYusef-Zadeh et al.1 (l2010h have recently noted an anticorrelation 
in NIR/X-ray emission with sub-mm/mm emission, as well in 
previously published data, extending the number of observed ex- 
amples to six. The same examples were in most cases originally 
interpreted as delayed sub-mm flares. This work suggested the 
anticorrelation could be due to an eclipse of the steady-state ac- 
cretion flow by the flaring region. 

In the following section, we present different modeling at- 
tempts to explain this interesting behavior observed in Sgr A* 
emission. 



4. SED models for an NIR/X-ray flare with a sub-mm 
dip 

VLBI measurements show that the size of Sgr A* de c reases 
from radio to the sub-mm (e.g. iBower et all 120041 120061 
iDoeleman et all 120081) . The sub-mm emission is universally 
produced in models by hot electrons of 6 e ~ 10 in the in- 
nermost regions of the accretion flow (close to the last sta- 
ble orbit), where the electron density is typically ~ 10 7 cm -3 



20031). In the GRMHD sim ulations of iDexter etall 



Moscibrodzka et al ] d2009l) the mm/sub -mm emis 



and magne t ic field strength is B~30G (e.g. lMarkoff et all 20011: 
lYuan et all 
(1201 Oh and 

sion region is situated near the midplane of the disk in the inner 
few Schwarzschild radii of the accretion flow. Though it is not 
at all clear where the NIR flares themselves originate, a dip in 
mm/sub-mm emission simultaneous to an NIR flare implies a 
causal connection between the two. In what follows, we discuss 
a number of possibilities for the relationship between the NIR 
and sub-mm emission region. We consider that the sub-mm dip 
might occur due to: 

- a shift in frequency of the sub-mm bump, with the sub-mm 
emission and NIR emission arising from low and high fre- 
quency sides of the sub-mm peak, respectively; 

- a decrease in the total number of sub-mm emitting electrons, 
through e.g. heating/acceleration (to produce the NIR flare) 
accompanied perhaps by cooling/escape of some fraction of 
the electrons; 

- a decrease in the sub-mm-e mitting electrons emissiv ity via a 
decrease in magnetic field (IDodds-Eden et al.ll2010l) ; 

- occultat ion of the quiescent com ponent by the flaring com- 
ponent dYusef-Zadeh et alll2O10h . 

We do not consider the expanding blob model here, for the 
simple reason that it does not predict that the highest NIR flux 
should be coincident with the lowest sub-mm flux. Furthermore, 
expanding blob models, which are presented elsewhere in the 
literature, have the following problems: (i) some observations 
appear to require a redder spectral index in the NIR than allowed 
by o ur lower limi t of f3 jj > -0.5 (see e.g. models A and B of Table 
4 in lEckart et all 120081) : the spectral index is given there as a 
where S v oc v~ a , so a — 1-/3, (ii) Moreover, they require self- 
absorbed flare emission in the sub-mm, which can only occur 
for the electron densities and/or magnetic fields that are much 
hig her than the typical value s in the accretion flow around Sgr 
A* (lDodds-Eden etalll201fj) . 
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Fig. 5. Flare model for Sgr A* from iFalcke & MarkofU d2000t) 
is plotted in blue and compared to the data from July 20 (purple 
square for the NIR average flux during the flare time interval and 
purple arrow for the MIR upper limit; note that the MIR upper 
limit from the July 22 was lower). The red stars present the peak 
flux of the July 20 flare, together with the simultaneous sub- 
mm flux. Blue circles show the median fluxes in a non-flaring 
state. The radio-sub-mm SED is shown in black (with squares: 
the compilation of Genzel 2010, diamonds: from Shcherbakov 
2011. The typical radio-sub-mm variability amplitudes and a 
com pilation of MIR u pper l imits are shown in gray and taken 
from lMelia&Falckd feoOll) . A flare due to a 30% increase in 
B in the jet model is also plotted as a red line. An increase in 
density resulting from a temporary increase i n the accretion rate, 
M, has the same result (iMarkoff et all 1200 ll) . 



4.1. The sub-mm bump 

A natural origin for an anticorrelation between NIR and sub- 
mm flux is that both high and low frequency parts arise from the 
same component: the sub-mm bump. One could imagine that if 
the peak of the sub-mm bump were to vary in frequency, e.g. due 
to underlying variations in density (th rough the accretion rate M 
or the magnetic turbulent as shown in lDexter et all |2010|) . while 
essentially preserving the peak sub-mm bump flux, one would 
observe an increase in the NIR flux as the sub-mm flux decreases 
(and vice versa). 

Such a setup is, however difficult to realize in practice. While 
an increase in density or magnetic field will increase the fre- 
quency at which the flux becomes self-absorbed, it also has the 
effect of increasing the flux at both sub-mm and NIR frequen- 
cies. This is illustrate d in FigUl which shows t he quiescent rel- 
ativistic jet model of IFalcke & MarkoffJ d2000h with a 30% in- 
crease in the magneti c field and an increas e in density produced 
by an increase in M dMarkoff et all l200ll) . This model was not 
created to fit our flare, but is shown more as a demonstration of 
the general behavior upon increasing B or n e in the steady-state 
models. There are two important aspects that are in conflict with 
our observations and that disfavor this model: (i) the sub-mm 
flux, as for the other frequencies, increases with the flare, there 
is no sub-mm dip; (ii) an exponential cutoff in the flare emission 
in the NIR is out of the question because of the MIR upper limit: 
nonthermal particles with at the most p = 4.6 are required. This 
reinforces that accelerat ion of nonthermal e l ectrons is needed 
to produce an NIR flare dMarkoff et all l200lHYuan et al.L l2003t 
iDodds-Eden et al.Ll2009h . 



4.2. Loss of electrons 



Another possibility for a decrease in the sub-mm emission is 
that a portion of the sub-mm-electrons were accelerated to NIR- 
emitting energies and/or were lost in an outflow. To determine 
how many electrons we need for the NIR flare, compared to how 
many we need to lose for the sub-mm dip, we create toy syn- 
chrotron models to model both the NIR flare emission as well as 
for the amount of "missing" sub-mm em ission. For the details of 
the models see Dodds-Ed en et al.l We model the missing 

emission with a thermal electron population of electron temper- 
ature 6 e = 10, magnetic field strength B = 30G, and electron 
density n e in a region of size R. We then model the flare with 
the same parameters, but with an accelerated electron distribu- 
tion, a powerlaw N(y) ~ y~ p of y m ,„ = 10 and p = 2.2 (since 
the numbers of electrons decrease so rapidly with y, the value of 
y cutoff does not have any effect on the electron density, as long 
as it is significantly higher than y m - ln and high enough to produce 
NIR emission). This choice of particle index both matches the 
published spectral indices for bright emissio n of B = 0.4 + 0.2 
(Hornstein et al 1 120071 iGflleTsen et all I2006T) . satisfies our MIR 
upper limits, and additionally ensures that the flare itself does 
not produce too much sub-mm emission (that would result in a 
simultaneous flare rather than a simultaneous dip). 

However, we find that the number of electrons required to 
produce the "missing sub-mm emission" in these models are 
greater than the number required to produce the flare, if the mag- 
netic field is the same before and after acceleration. Only 30% 
of these sub-mm-emitting electrons are required to produce the 
NIR flare under the conservative assumption y, m „ = Qe- This im- 
plies that the missing sub-mm emission cannot be explained by 
the acceleration of sub-mm-emitting electrons alone: some frac- 
tion must be completely lost to the system in the process. 

Given our assumption that electrons are accelerated out of 
a thermal population, this result is surprisingly robust to the 
choice of parameters: i.e. 6e, B, for the unaccelerated popula- 
tion of electrons and y m ,„, p for the accelerated population. This 
is because we are essentially only looking at a relative change 
in the electron distribution. Both the missing sub-mm emission 
and the flare emission, for example, are subject to the same B by 
assumption, so changing B essentially only requires a rescaling 
of the density (which is, however, the same for both unaccel- 
erated and accelerated distributions, so a similar fraction must 
still escape). The fraction of escaping electrons is more sensi- 
tive to y m j n of the accelerated population. A lower y„„-„ however 
is unrealistic under the assumption of electron acceleration: our 
condition of y ;)1 ,„ = Qe is already quite conservative since only 
«7% of electrons actually have energies lower than Qe in a rel- 
ativistic Maxwellian distribution (the typical energy is, in fact, 
y = 29e); in other words, electrons with y < Be might as well 
be treated as escaped for the purposes of this paper. A higher 
ymin > Be only decreases the number of electrons required for 
the flare (since the density depends sensitively on y„„„), thereby 
increasing the fraction that mus t escape. Since it is quite likely 
that y mm > 6e dDing et all 1201 Oh . we might expect a much more 
significant fraction of electrons to escape. These electrons could 
be carried away in a very fast outflow, for example, due to their 
small number producing very little emission as they rapidly ex- 
pand and cool (adiabatic expansion IDodds-Eden et ail, 1201 Oh or 
may be accreted onto the supermassive black hole. 
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4.3. Decrease in magnetic field 

In lDodds-Eden et all d2010h . it was suggested that a decrease in 
magnetic field (for instance, that accompani es a magnetic recon- 
nection event or due to magnetic turbulence Dexter et al.l 1201(H) 
could cause a decrease in the longer wavelength emission. This 
could be more efficient at producing a dip in flux than loss of 
electrons since the synchrotron power emitted by an electron is 
proportional to B 2 , while it is only linearly proportional to the 
number of electrons. 

In our toy model, we now hold the total number of electrons 
constant (no escape) between the thermal (Of ~ 10) and accel- 
erated power-law models (y m ,-„ = 10, p = 2.2) described in the 
last section. The accelerated power-law model, with the same 
number of electrons as the thermal model and the same mag- 
netic field strength, overproduces the NIR emission. Instead of 
decreasing the number of electrons, we find that if we decrease 
the magnetic field from B = 30G to B — 14G, the observations 
are well-matched. This flare mode l is sho wn in Fig J6] compared 
and combined with the I Yuan et al.l (120031) model whose thermal 
component was scaled down to match the SED in the sub-mm 

Sart of the section (in particular our blue "no flare" data point) 
. As in the case of electron loss, y„„„ = Of is quite conserva- 
tive; if the true y mm is higher we would have a higher density 
of high-energy accelerated electrons, requiring an even greater 
decrease in the magnetic field to prevent overproduction of NIR 
emission. However, the magnetic field in a realistic magnetic re- 
connection event probably als o decreases over a larger region 
than the acceleration site itself dDodds-Eden et al.ll2010l) . which 
would produce a larger sub-mm dip, while allowing the freedom 
of a more moderate decrease in B. 

4.4. Eclipse 

In lYusef-Z adeh et al] (1201 Oh it was proposed that a decrease in 
sub-mm emission could occur via an eclipse of the steady-state 
emission by the temporary flaring blob. In this section we ex- 
plore the possibility that occultation is responsible for the sub- 
mm dips observed in this paper. In this model it is considered 
that the two regions are separate: i.e. accelerated flare electrons 
do not come from (or produce very little disturbance to) the qui- 
escent population, and the flaring region is in front of, not em- 
bedded in, the quiescent region in the line of sight. 

Fig ur j7] shows a plot of the equations of lYusef-Zadeh etail 
(1201 Oh for variation in sub-mm flux vs. size of the flaring region 
(Rf), for different quiescent region sizes (Rq). Anything larger 
than Rq > IRs encounters very little occultation. Evidently, both 
the 850/mi and NIR emission regions would have to be very 
small (Rq < 0.6Rs and Rf = 0.5-1.0 Rq) to produce a dip as 
large as 0.5 Jy in the sub-mm emission via occultation. 

A quiescent size of 0.5Rs i s much sm a ller th an would be ex- 
pected. Using the caveat from iMarscherl (119831) (radius = 0.9* 
FWHM ) on the FWHM size 3.7^ lS (37 microarcseconds) at 
1.3mm (iDoeleman et al.L 120081) an d extrapolating with the ra- 
dius oc A 13 wavelength dependence iBower et al.l (120061) . we ob- 
tain a radius at 870/im of 2.0 Rs- However, extrapolating the 
size downwards from 1.3 mm to 850 pm could be incorrect as 
the accretion flow might already become optically thin around 
1.3 mm. In this case the size would no longer change. 



5 W e note that the nonthermal component of the model of lYuan et alj 
with p = 3.5 lies close to the non-flaring A".s-band flux of Sgr 
A*, which is als o similar to the long-term median of the continuously 
variable NIR flux lDodds-Eden et alJooTTh . 
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Fig. 6. A model for an NIR flare with an accompanying sub- 
mm dip compared to the data from July 20 (see caption of 
Fig|5]l. Compared to the sub-mm-radio SED shown in black 
square/diamonds, it can be seen that the flux of the sub-mm dip 
(red star) is lower than the typical sub-mm flux. The dashed lines 
show a toy model for the flare, in red, as well as for the "missing" 
sub-mm flux, in blue. The blue solid line shows the quiescent 
model o f lYuan et al.l (120031) . and the red solid line represents the 
same model but with the "missing" sub-mm flux subtracted and 
the flare emission added. This particular figure shows the de- 
creasing magnetic field case: to prevent overproduction of syn- 
chrotron emission from the accelerated population the magnetic 
field decreases from B = 30G to B — 14G. The spectrum for the 
case of the loss of electrons is very similar. 



Moreover, the corresponding flare size would be so small 
that it would imply a very high density of nonthermal parti- 
cles. It is possible that relaxing the assumption o f equipartition 
in the eclipse model of lYusef-Zadeh etail ( 120101) might change 
these conclusions. Another possibility is that the electron spec- 
tral index might be steepe r than the assumed p — 2 of the model 
( lYusef-Zadeh et aUl2010b . Although a steeper electron spectrum 
can increase the optical depth at low frequencies (and thereby 
the absorption of quiescent emission), it also increases the sub- 
mm emission from the flare itself significantly (on the order of 
Jy; iDodds-Eden et all l2010h . which is likely to overwhelm any 
increased absorption of the quiescent sub-mm emission. 

Some radio ( 43 GHz) dips have been reported as well dur- 
ing an NIR flare (lYusef-Zadeh et al.Lp2010h . We note that, for the 
same modeling parameters we used, no significant dips in the ra- 
dio are expected for any Rf , given a 43 GHz (i =7mm) emitting 
region size of R Q = 11R S (IBower et all 120061) . Only a smaller 
(43 GHz) quiescent region size can produce observable radio 
dips due to occultation, for example for Rq = 4R$ (four times 
smaller than observed) radio dips of up to ~ 20% of the emission 
could be produced for flare region sizes of 1 - 2R$ . Yusef-Zadeh 
et al. argue that a 4x smaller quiescent region size was not neces- 
sarily in conflict with the size radio measurements as the source 
might be much smaller in semi-minor axis; however, this in turn 
imposes increasingly difficult constraints on both the geometri- 
cal shape of the flaring region and its alignment with respect to 
the quiescent region to retain a high degree of occultation (given 
that a covering fraction (Rf/Rq) 2 of at least 0.5 2 is required). 
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0.1 1 
R F [Rs] 

Fig. 7. Dips in sub-mm emission via occultation by a flaring re- 
gion of size R F , given a sub-mm flux of f 350GHz = 3 Jy, appro- 
priate for our observations of July 20. The different lines show 
the size of the dip that can be expected for different quiescent 
region sizes, Rq. 

5. Conclusion 

We have reported on a sub-mm, NIR, and MIR observation cam- 
paign of Sgr A* over five nights in July 2007. First, we de- 
tected NIR flare events of different amplitudes and structures. 
We also estimated the lowest upper limit ever derived for a flare 
at 8.59//m (22.4 mJy with = 1.6 + 0.5) and derived 

constraints on the spectral slopes between NIR and MIR. The 
analysis of the correlation between sub-mm and NIR variabil- 
ity shows that the highest NIR brightness levels are coincident 
with the lowest sub-mm levels. We examined several possible 
explanations for this observational fact. We disfavor variations 
in the sub-mm bump mainly because the sub-mm and NIR emis- 
sion from only one thermal component does not satisfy our MIR 
upper limit. Occ ultation of the flaring com ponent by the quies- 
cent component (lYusef-Zadeh et all , l2010l) seems unlikely given 
the small sizes of the emitting areas at play. We argue that the 
NIR/sub-mm anticorrelation could likely be due to either elec- 
tron loss (>70% for y m j„ — 6e ), a decrease in the magnetic field 
(from B = 30G to B — 14G), or a combination of both, involv- 
ing fewer escaping particles and a smaller decrease in magnetic 
field. It is worth noting that a decrease in magnetic field, as ex- 
pected during a magnetic reconnection event, can also explain 
the differen t durations of observed fl ares in the NIR and X-ray 
lightcurves dDodds-Eden et all I2.Q10I) . Regardless of the mecha- 
nism, we conclude that there is a link between a decrease in sub- 
mm emission and the acceleration of nonthermal electrons (or 
thermal electrons with 6>£ >2000). While more work is needed 
to pin down the nature of this link precisely, magnetic activity, 
fast outflows, and their interaction with the black hole are very 
likely to be pieces of the puzzle. 
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